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ABSTRACT 

We used the GEOS database to study the Blazhko effect of galactic RRab stars. The database 
is continuously enriched by maxima supplied by amateur astronomers and by a dedicated survey by 
means of the two TAROT robotic telescopes. The same value of the Blazhko period is observed at 
different values of the pulsation periods and different values of the Blazhko periods are observed at 
the same value of the pulsation period. There are clues suggesting that the Blazhko effect is changing 
from one cycle to the next. The secular changes in the pulsation and Blazhko periods of Z CVn 
are anticorrelated. The diagrams of magnitudes against phases of the maxima clearly show that the 
light curves of Blazhko variables can be explained as modulated signals, both in amplitude and in 
frequency. The closed curves describing the Blazhko cycles in such diagrams have different shapes, 
reflecting the phase shifts between the epochs of the brightest maximum and the maximum O-C. Our 
sample shows that both clockwise and anticlockwise directions are possible for similar shapes. The 
improved observational knowledge of the Blazhko effect, in addition to some peculiarities of the light 
curves, have still to be explained by a satisfactory physical mechanism. 

Subject headings: Astronomical data bases: miscellaneous - Stars: evolution - Stars: horizontal-branch 
- Stars: variables: RR Lyrae 



1. INTRODUCTION 

In the early XX th century S.N. Blazhko pointed out 
cycle-to-cycle variations in the light curves of Var 87.1906 
Draconis=RW Dra, di scovered by L.P . Tserasskaya at 
Moscow Observatory (jTsesevichl 119691 ). The periodic 
changes in the times of maximum brightness were fitted 
by adding a sinusoidal term to the linear ephemeris, as 
happened a few years later when t he same effect was ob- 
served in the case of RR Lyr itself (lShaplevlll916h . After 
being discovered in several other RR Lyr stars, the vari- 
ations of the times of maximum brightness (T max ) and 
of the magnitudes at maximum (V max ) were defined as 
the Blazhko effect since the list of maxim a recorded from 
July 15 to August 25, 1906 on RW Dra (lBlazhkdll907l) 
constituted the first evidence of it. iKlepikoval (|1956f) re- 
ported on one of the first observational projects to study 
the Blazhko effect in a systematic way. 
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The characteri stics of the Blazhko stars disco vered in the 
XX th century ()Szeidll97d[T^lSmitI]|[T995] ) were deter- 
mined by searching for periodicities in the V max and/or 
O-C (observed minus calculated T max ) values and/or by 
performing a frequency analysis of the photometric time- 
series. In the latter case the periodic variations of the 
shape of the light curve cause a multiplet centered on 
the pulsational frequency and with a separation equal to 
the Blazhko frequency. 

The number of Blazhko stars did not increase signif- 
icantly until the advent of CCD techniques. Stars 
much fainter than the threshold for photolectric pho- 
tometry could be discovered in large-s cale surveys 
as M ACHO (MAssive Compact Halo; lAlcock et all 
120031) a nd OGLE (Optical Gr avitational Lensing Exper- 
iment; iSoszinski et al.l 120101 and references therein). 
iMoskalik fc Porettil (|2003t) suggested that -23% of the 
RRab stars in t he ga lactic bulge shows a Blazhko effect. 
ICollinge et al.l (|2006l ) substantially confirmed this inci- 
dence (27.6%) by performing an analysis on a larger sam- 
ple and a longer baseline. A much larger incidence was 
found in the Konkoly Bl azhko Surveys (KBS s): 14 out of 
30 (47%) in the KBS I (lJurcsik et al.ll2009bD an d 45 out 
of 105 (43%) in the KBS II (|S6dor et al.ll20ll . This is 
probably due to the very small modulations which could 
be detected thanks to the multicolour CCD photome- 
try performed with an automatic 60-cm telescope. The 
Konkoly surveys also provide information on the changes 
in the physical parameters during the Blazhko cycle. 
The observational scenario changed with the continu- 
ous monitorin g of Blazhko RR var i ables with the satel- 
lites CoRoT dchadid et al.l [2011 iPoretti et aO [20T0I: 
Guggenberger et al.l 120111) and Kevler (|Kolenberg et al.l 



2010; Sza bo et al.l 12010). The lack of any relevant alias 



structure in the spectral windows of the space-based 
timeseries allowed us to obtain new results, as evidence 
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of large cycle-to-cyc l e var iations of the Blazhko effect 
(|Guggenberger etall 120111) of the excitation of addi- 
tional modes ()Poretti et all 120101). an d of the period dou- 
bling bifurcation (jSzabo et al 

Despite these new observational results, the physical ex- 
planation of the Blazhko effect still remain controversial. 
Very differen t mechanisms were invoked: th e oblique pul- 
sator model (lKurtzHl982t iShibaha shi 2000 ]), the resonant 
pulsator model (jDziembowski fc Mize rski 200 41) , and th e 
action of a turbulent convective dynamo (|Stothersi r2006). 
None of these models is able to p rovide a close matching 
between theory and observations ([Ko vacs 2009) and then 
it seems that we are still far from a satisfactory theoret- 
ical explanation. 

The GEOS (Groupe Europeen d 'Observations Stellaires) 
RR Lyr Databasc(3 aims at collecting the times of max- 
imum brightness of RR Lyr stars, exclusively galac- 
tic RRab (fundamental radial pulsators) and RRc (first 
overtone radial pulsators). The project was started in 
2000 by archiving recent and historical publications of 
T max in the literature. It now contains 61000 maxima of 
3600 stars and it is regularly fed by the valuable inputs 
from amateur astronomers (BAV, AAVSO, GEOS ...). 
We already used the GEOS database to obtain rele- 
vant results in the study of the evoluti onary changes of 
the p ulsational periods of RRab stars (|Le Borgne et al.l 
|2007| ). The same data can be used to study period vari- 
ations for a wide sample of galactic RR Lyr stars on a 
much shorter timescale, in order to bring out new tiles 
necessary to compose the puzzle of the Blazhko effect. 

2. THE GEOS RR LYR SURVEY WITH TAROT 
TELESCOPES 

The main contribution to the GEOS RR Lyr database 
comes from the GEOS RR Lyr Survey and particu- 
larly from the observations of the two robotic telescopes 
TAROT (T elescope a Action Rapid e pour les Objets 
Transitoires; iKlotz et al.l I2008L [2009) located in France 
(Calern Observatory) and in Chile (La Silla Observa- 
tory). These 25-cm automated telescopes are dedicated 
to the capture of the optical afterglows of Gamma Ray 
Bursts. Several programs are scheduled in the idle time, 
as the one on galactic RR Lyr stars. Indeed, RR Lyr 
stars are very suitable targets for small robotic telescopes 
due to their short-period and large-amplitude variability, 
which allow the observation of several T max and V max on 
a single night. TAROT telescopes have produced about 
840 max/year since the beginning of the observations in 
2004. This is about half the total production of 1640 
max/year from ground-based observations on the same 
period. This performance appears more relevant when 
compared to an average of about 500 max/year over the 
110 years elapsed since the discovery of RR Lyr stars. 
RRab stars have periods in the 0.37-0.65 d range and 
we could collect from each TAROT site no more than 
1 maximum/night for each given star. A few equato- 
rial stars could be observed both from Calern and La 
Silla, but observing two consecutive maxima is an excep- 
tional event due to weather differences. Since the main 
goal is the monitoring of as many RR Lyr stars as possi- 
ble, the scheduling software tries to secure points for all 

12 The GEOS database is freely accessible on the internet at the 
address http : //rr-lyr . ast . obs-mip . f r/dbrr/ 
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Fig. 1. — Spectral windows of three RR Lyr stars having different 
periods. 



the program stars and hence gives a lesser priority to a 
star that has already been observed on a previous night. 
The telescopes are continuously moved from one star to 
the next during the night and the measurements consist 
of two consecutive 30-s exposures taken every 10 min. 
They are performed around the predicted T max and the 
ephemerides are continuously updated to avoid to miss 
the real maximum. The mean error bar of the observed 
T max is about 0.003 d (4.3 min). We obtain an average 
of 50 measurements in 4.8 hours for each maximum. 
The resulting spectral window clearly shows alias peaks 
at ±1 y" 1 = ±0.0027 d" 1 (Fig.[TJ inserted boxes), due to 
the pointing limitations when stars are too close to the 
Sun, and different alias structures appear as a function 
of the period. The maxima of the shorter-period RR 
Lyr stars could be observed very frequently and their 
spectral windows are quite clear, as shown by that of 
AH Cam (Fig. Q] top panel). Stars with a period close 
to 0.50 d could not be observed for a long time when 
the maxima occur around sunset or sunrise and the gaps 
in the timeseries produced alias peaks in their spectral 
window (RZ Lyr, middle panel). If the period is around 
or longer than 0.6 d the next observable maximum af- 
ter a positive determination occurs a few days later (at 
least two). Therefore, the pseudo-Nyquist frequency 
is around 0.25 d _1 , as shown by the case of ST Boo 
(lower panel). This extreme case is still largely satisfac- 
tory for our purposes, since the shortest k nown Blazhko 
perio d (Pb) is that of SS Cnc, 5.309 d (|Jurcsik et al.l 
2006), corresponding to Jb =0.188 d _1 . Therefore, we 
performed our frequency analyses in the interval 0.0- 
0.2 d _1 b y means of the iterative sine-wave least-squares 
method (|Vaniceld[T97l . Due to the scheduling proce- 
dure of the TAROT telescopes, photometric measure- 
ments are concentrated around the maxima. The rising 
branch of the light curves are well covered including the 
shoulder in some stars currently interpreted as the ef- 
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feet of a h ypersonic shock-wav e propagating in the at- 
mosphere (|Chadid et al.l [20081 ) . Our analysis detected 
such a feature in the light curves of the non-Blazhko 
stars BC Dra and AR Per (Fig. [2]). BC Dra is a long- 
period RRab star (P=0.72 d), AR Per a short-period 
one (P=0.42 d). The almost complete phase coverage 
of the BC Dra curve allowed us to determine a pul- 
sational period of 0.71958083±0.00000015 d. The am- 
plitudes of the harmonics decrease up to 7/ and then 
they show small fluctuations before the final decline at 
12/. The residual rms is 0.033 mag and the residuals still 
shows the sudden change in the light curve occurring on 
the rising branch (Fig. [21 lower panel). The frequency 
analysis of the TAROT timeseries and the 105 observed 
maxima did not detect any signature of the Blazhko ef- 
fect. We applied the same procedure to the AR Per 
data, thus obtaining P=0.42555060±0.00000010 d in the 
2004-2011 interval. This accurate determination is par- 
ticularly useful s ince the period of AR Per is subject to 
secular changes (|Le Borgne et al.l 120071 ) . As a drawback 
of the TAROT scheduling procedure, the phase cover- 
age of the light curves of many program stars is often 
incomplete, with large gaps on the descending branch. 

3. THE TEST BENCH: STARS WITH A KNOWN 
BLAZHKO EFFECT 

Due to pointing accuracy and small seasonal effects 
affecting it, the reduction software cou ld select differ- 
ent s tars for the ensemble photometry (jDamerdji et al.l 
2007). This translates into small systematic differences 
that in some cases could prevent the frequency analy- 
sis of the V max values. The T max values are free from 
this effect, since it is independent from the choice of the 
reference stars. However, the O-C values depend on the 
stability of the main pulsational period, which could also 
vary on a time interval as short as seven years. A linear 
or parabolic ephemeris co uld locally fit the T max values 
in an unsatisfactory way (|Le Borgne et al.l I2007T) . leav- 
ing some residual peaks at very low frequencies. In such 
cases we corrected this secular effect by calculating a lin- 
ear ephemeris in the interval covered by the T max values. 
Finally, we considered both O-C and V max values to de- 
tect and characterize the Blazhko effect, giving full con- 
fidence to the results only if they were in agreement with 
each other. With the exception of just a couple of cases 
which will benefit from a longer time baseline covered 
by accurate maxima determinations, this double-check 
ruled out the possibility that the periodicity observed 
in the O-Cs was caused by a light-time effect (i.e., a RR 
Lyr variable in a wide binary), since V max values are not 
expected to show any appreciable variation due to the 
orbital motion. We have a couple of cases where a long- 
period, light-time effect is plausible: TAROT telescopes 
are monitoring them for a decisive confirmation. First, 
we analyzed the stars for which the Blazhko effect was 
already known. Figure [3] gives an overview of the tools 
applied to some stars that will be discussed later: power 
spectra of O-C and V max values (for RR Gem and AH 
Cam); the folded curves of O-C and V max values (for RR 
Gem and SZ Hya); the folded curves of original measure- 
ments (for SZ Hya); and the folded curves of residuals 
(for AH Cam). As shown in the upper part of Table Q] 
our independent analysis confirmed most of the previ- 
ously known results. For instance, the very small am- 
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plitude Blazhko effect of RR Gem ([Jurcsik et al.ll2005af ) 
was correctly identified and characterized (Fig. O top 
panels). This not only strengthened the confidence in 
our approach, but also allowed us to build the first O-C 
vs. V max plots for these stars. Without entering in a 
detailed discussion of the well-established cases, we em- 
phasize here some methodological and scientific items. 

3.1. The light curves and the periods of BD Dra and 
AH Cam 
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Fig. 4. — Blazhko effect of BD Dra. Left panels: power spectra 
of O-C and V max values {top) and curves of the O-Cs and V m ax 
values (bottom). Right panels: TAROT CCD measurements (top) 
and residuals after the 24-frequency fit (bottom). 



Fig. 5.— Blazhko effect of XY Eri. Top panel: TAROT CCD 
measurements folded with the pulsational period Po=0. 554261 d. 
Bottom panels: power spectra (right) and folded curves (left) of 
O-C and Vmax values. 



The dispersion in the residual plot was also observed 
on the rising branch of the Blazhko stars BD Dra and AH 
Cam. Regarding BD Dra, the analysis of 44 old photo- 
graphic and visual maxima did not show any significant 
periodicity. On the other hand, the frequency analysis of 
the 136 CCD maxima collected from 2005 to 2011 imme- 
diately pointed up a very clear periodicity of 24.11 d in 
both the O-C and V max values (left top panel in Fig.[4|. 
The O-C and V max values folded with the Blazhko period 
(P B =24.107±0.001 d) show a full amplitude of 0.043 d 
and 0.40 mag, respectively (left bottom panel in Fig.|4|. 
The Blazhko effect is clearly noticeable in the folded light 
curve (P o =0.58900964± 0.00000010 d; right top panel in 
Fig. [4]). Very sharp triplets were also detected around 
fo and harmonics up to 7/o- It was not possible to in- 
vestigate the existence of a quintuplet structure, since 
such peaks do not stand out clearly over the noise level 
of 0.003 mag (e.g., the component /o + 2 fs has an ampli- 
tude of 0.009 mag only). The plot of the residuals after 
a fit with 24 terms (n/o ± /s, with n=l,..., 8) shows a 
large spread in the values at the phases corresponding 
to such part of the light curve (right bottom panel in 
Fig. ED. 

The TAROT maxima collected on AH Cam in the recent 
years clearly indicate that the true Blazhko frequency 
is / B =0.092341 d" 1 (Fig. El middle panels), i.e., the 
alias at +1 y _1 of the previous value (|Le Borgne et al.l 
2007). This improvement in the analysis is due to the 
long coverage ensured by the automatic telescope during 
each year, able to strongly damp the aliases at ±1 y -1 . 
The Blazhko frequency was also found in the analysis 
of the TAROT original measurements as triplets cen- 
tered around /b=2.71211 d _1 and harmonics up to 5/o- 
The residuals show a large scatter in the rising branch 
(right-middle panel in Fig. [3]) , although not as large as 
in the case of BD Dra. The periods obtained from the 



frequency analysis are P=0.36871596±0.00000006 d and 
P B = 10.8289±0.0002 d. 

3.2. The O-C and Vm ax amplitudes of the Blazhko 

effect 

The star with the most spectacular Blazhko effect is 
XY Eri (Fig. [5]). The V max values vary in a range of 
0.75 mag, more than half the full amplitude of the peak- 
to-peak amplitude, and the O-Cs in a range of 91 min. 
The Blazhko effect of IK Hya is very similar, since the O- 
Cs have the same range, but the V max amplitude is prob- 
ably around 0.5 mag. Since we only observed 9 maxima, 
IK Hya deserves further investigation in the future. Most 
of the stars shows a V max amplitude between 0.20 and 
0.40 mag and only RS Boo, SS Oct, and RR Gem have 
amplitudes smaller than 0.20 mag. These large variations 
explain why the Blazhko effect was sometimes referred to 
an amplit ude modul ation (see also the introductory re- 
marks in Szeidl fl98ct ) . Note how the O-C values show 
variations only in a narrow phase of the Blazhko period 
of SZ Hya (Fig. [3l bottom panels). SS CVn shows a sim- 
ilar variation of the O-Cs, but its Blazhko period is 3.6 
times longer. 

3.3. Long Blazhko periods and long-term modulations 

Blazhko periods longer than 100 d are less common 
than shorter ones. In our sample we have the cases of 
RS Boo, ST Boo, RZ Lyr, RX Cet, and UV Oct. The 
longest Blazhko period is that of RS Boo (533 d). The 
modern CCD observations provide a spectacular confir- 
mation of the second longest Blazhko period, that of ST 
Boo (Fig. [6]), whose features were never presented before. 
The analysis of the TAROT CCD timeseries of RZ Lyr 
supplies P=0.511241±0.000005 d a nd P B =120.19 d, thus 
confirming the trends reported by Uurcsik et al.l (|2012[ ) 
on a longer time baseline. In the case of RX Cet, the 
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Fig. 6. — Blazhko effect of RW Dra and ST Boo. Power spectra 
and folded curves of the O-C and V max values of RW Dra (left 
panels) and ST Boo (right panels). 

15 maxima spread over 7 y did not aiiow a precise de- 
termination of the Blazhko period, characterized by the 
second largest amplitude (0.061 d) of the O-Cs after that 
of XY Eri. We detected a tentative period of 273 d su- 
perimposed on a slight secular variation of the period. 
On the other hand, UV Oct shows a low O-C amplitude 
of only 0.012 d. 

Long, additional modulations were also claimed as sec- 
ondary periodicities superimposed on sho rt Blazhko pe- 
riods. For instance, iKholopov et al.l (|1985l) report on a 
composite Blazhko effect for RW Dra. The O-C values in 
the GEOS database clearly indicate the Blazhko period 
of 41.4 d (Fig. |6j). The power spectrum and the folded 
light curve of the residual O-C values are not support- 
ing the action of other strictly periodic cycles. On the 
other hand, the V max values show some scatter, thus sug- 
gesting changes in the amplitude from one cycle to next 
(Fig. HJ . XZ Dra is perhaps the extreme case of changes 
in the Blazhko cycle. We did not detect any reliable 
periodicity around 75 d in the recent maxima collected 
in the GEOS database. The dampening of the Blazhko 
variation has been noticed for mo re than 50 years in the 
XX th century (jJurcsik et al.ll2002l ). Since still more rapid 
variations are clearly detected in t he continuous space 
times eries (e.g., CoRoT 105288363, Guggenb erger et al.1 
1201 if ), we argue that the composite Blazhko effect and, 
more generally, the long-term variations could be a spu- 
rious effect of such cycle-to-cycle changes, irregularly or 
poorly sampled from ground-based monitoring. 

3.4. Change in the Blazhko period with changes in the 
pulsational period: Z CVn 

We detected the Blazhko effect in the CCD data of 
Z CVn, together with a well-defined secular change 
of the pulsation period. The pulsational period was 
0.653942±0.000007 d around JD 2439000. The more re- 
cent CCD data were obtained in a time interval where 



Fig. 7. — Changes in the Blazhko effect of Z CVn. Top panel, 
left: the plot of the O-Cs with respect to a linear ephemeris shows 
a secular change in the period. Top panel, right: folded TAROT 
light curve showing the Blazhko effect. Bottom panels, left: the 
power spectra of the O-C and V ma x values clearly show the peak 
a t /b— 0.044 d — 1 . Bottom panels, right: O-C and V ma x 

folded 

with the Blazhko period corresponding to fg. 

the pulsational period has undergone a sudden change 
to 0.6538907±0.0000003 d after JD 2454200 (Fig. [2 top 
panel). The Blazhko period around JD 2439000 was 
22.75 d (| K an vol [I 9 6 61 ). This value is confirmed by our 
re-analysis of the T max values, even if we have to note 
that the power spectrum shows a broad peak. When 
analysing the CCD data we immediately noticed that 
the amplitude of the O-C variation was 0.021 d after 
JD 2454200, while it was 0.033 d in Kanyo's data. This 
is the first remarkable change in the Blazhko modulation. 
More important, the frequency analysis of the two sets 
of O-C and V max , values supplied a Blazhko period of 
22.95±0.03 d after JD 2454200 (Fig. [7] left-bottom pan- 
els). The amplitude of the O-C values, corrected for the 
period change, decreased to 0.021 d, while the amplitude 
of the CCD V max values remained constant at 0.32 mag 
(Fig. [3 right-bottom panels). 

Our analysis shows that the pulsational period affects 
the Blazhko one. In particular, the decrease of the pulsa- 
tional period resulted in a longer Blazhko period, i.e., the 
two changes are anticorrelated. T he same behaviour wa s 
observed in the cases of XZCyg (ILaCluvze et al.l[ 2004). 
RV U Ma dHurta et al.l I2008D. R W Dra (iFi rmanvukl 



I1978D . and RZ Lyr (jJurcsik et al.ll2012fl . However, this 
does not seem a strict rule, since correlated changes were 
observed in the case of XZ Dra (jJurcsik et al.ll2002l ). 

4. NEW BLAZHKO STARS 

Corroborated by the results obtained on stars with a 
known Blazhko effect, we applied the same tools to stars 
poorly observed in the past. TAROT photometry was 
able to detect the signatures of the Blazhko effect in the 
timeseries of the stars listed in the lower part of Table [TJ 
The Blazhko frequencies are determined without ambi- 
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Fig. 8. — Stars with the brighest V max (ordinate) corresponding to the maximum positive O-C (abscissa). The closed curves are run 
anticlockwise. 
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guity in the cases of SS CVn, TY Aps, BI Cen, AF 
Vel, and SV Vol. The Blazhko effect of ST Vir is 
very small, similar to that of RR Gem, and its de- 
tection is another good example of the effectiveness 
of the TAROT survey. In the cas e of B I Cen, the 
uncertain value reported by iSmithl (| 19951 ) (70 d) is 
slightly different from what we determined. The light 
curve was well covered by TAROT observations and 
we could determine P=0.45319439±0.00000007 d and 
Pb=79.75±0.02 d from the frequency analysis of the 
timeseries. BI Cen shows very sharp O-C and V max vari- 
ations and the Blazhko variability looks like the most 
regular in our sample. 

The fact that most of the timeseries are single-site in- 
troduced some uncertainties on the determination of the 
true Blazhko frequencies, which could be actually the 
aliases at ±1 y" 1 or ±2 y^ 1 of the highest peak in the 
power spectra. This ambiguity results in a real period 
shorter or longer by about 10 d for a Blahzko period of 
40 d (BS Aps and UU Hya) and 4 d for a period of 23 d 
(U Cae). The situation is quite different 
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Fig. 9. — Same as Fig. [8] but the closed curves are run in the 
clockwise direction. 

for the stars showing a long Blazhko period, i.e., RX Col 
and RY Oct. At frequencies of 0.008 d _1 (correspond- 
ing to a period of 130 d, as that of RX Col) or lower 
(RY Oct) the alias at —1 y _1 points to an almost twice 
longer period. The extreme case is that of RU Cet. The 
O-C's, the V max , and the folded light curve suggest the 
possibility of a long period, probably much longer than 
100 d, but the strong aliasing effect prevented us from 
proposing a reliable value. 

In addition to those listed in Table [TJ other RRab stars in 
the GEOS database show cycle-to-cycle variation which 
could be ascribed to a Blazhko effect, but we could not 
propose a reliable value for the period due to the strong 
alias uncertainty. The scientific discussion of the Blazhko 
effect could be sharpened if new values for T max are col- 
lected on these stars. 

5. DISCUSSION 

5.1. Varieties in the 'phenomenology of 'the Blazhko 

effect 

With respect to the frequency analysis of ground- or 
space-based photometric timeseries, the GEOS database 
enriched with the maxima observed with the TAROT 
telescopes offers the possibility of investigating the re- 
lation between amplitude and p hase changes, descr ibed 
by the O-C vs. V max diagram. IBenko et al.l (|2011| ) re- 
cently interpreted the Blazhko light curves as modulated 



signals and provided theoretical predictions of these dia- 
grams (see their Figs. 17 and 18). 

To exploit this possibility, we used the O-C and V max 
folded with the Blazhko period (see bottom left panels in 
Fig. HI for an example). We calculated the least-squares 
fits of these folded data, thus obtaining an O-C, V max 
couple at any given Blazhko phase. We drew the closed 
curve that describes the shape of the Blazhko effect by 
connecting the consecutive couples. A quick analysis im- 
mediately reveals how such closed curves can assume dif- 
ferent shapes and can be run in opposite directions. We 
have 20 cases of anticlockwise directions and 9 cases of 
clockwise ones. The most numerous category is formed 
by Blazhko cycles where the O-C and V max curves are 
roughly in phase and the closed curves are run in the 
anticlockwise directions (13 cases, Fig. [5]). The A</> shifts 
range from -0.06 P B (RS Boo) to -0.26 P B (TT Cnc). 
The curves of Z CVn (A0 = -0.79 P B ) and UU Hya 
(—0.83 P B ) have the same orientation, but run clockwise 
(2 cases, Fig. [9]). 

The O-C and V max curves are roughly in anti-phase 
and the closed curves are run anticlockwise for six stars 
(Fig. [TO]). The shifts range from -0.30 P B (RX Col) 
to —0.42 Pb (UV Oct). For this orientation of the closed 
curve the number of stars running clockwise is nearly the 
same (seven cases, Fig. [Til and the phase shifts ranges 
from -0.52 P B (DD Hya) to -0.66 P B (RZ Lyr). 
There are close similarities with the curves predicted by 
the a mplitude and frequency modulations (Ben kd et al.l 
120111) . In particular, the extreme cases of the Blazhko di- 
agrams shown by BD Dra, SZ Hya (Fig. [HJ), and VW Dor 
(Fig. ITTj) confirm the complicated shapes predicted in 
case of simultaneous non-sinusoidal amplitude and phase 
modulations. They are characterized by strongly elon- 
gated shapes, with some cuspids and loops, due to the 
peculiar plots of the O-C and/or V max values. In the 
case of SZ Hya, (Fig. [3] and Fig© the A0 shifts reflect 
this peculiarity and is different from the others of the 
same group of stars. 

These subdivisions reflect the mat hematical difference s 
in the signal modulations stressed by Benkd ct"ail (|2011|) . 
The phase shifts A</> (Table [IJ can be tranf ormed in the 
phase parameter (f> m (Eq. 45 and Fig. 17 in [Benkd et al.l 
1201 ID by the relation 4> m = A<j) + k 180, where k = 1 
for the anticlockwise direction, fc = — 1 for th e clock wise 
direction. It is noteworthy that IBenko et al.l f|201 ID ob- 
tained their diagrams by a mathematical representation 
of the light curve, while ours were built by using the ob- 
served T max and V max values, without any assumption 
on the form of the light curve. We can argue that our 
observational results are a rather solid confirmation that 
the Blazhko effect can be described by amplitude and 
frequency modulations, though the physical reasons of 
these modulations still remain unknown. 

5.2. Distribution of the Blazhko periods 

The characteristics of the all-sky galactic Blazhko RR 
Lyr stars could be comp ared with those obta ined in the 
galactic bulge by OGLE (Collin ge et al.ll2006[) and in the 
Magellanic Clouds by MACHO ( jAlcock et al.ll2003l ). The 
results of the OGLE and MACHO surveys were obtained 
by performing a frequency analysis of the photometric 
timeseries and the Blazhko periods were measured from 
the separation of the multiplets in the amplitude spec- 
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TABLE 1 

RRab stars showing Blazhko effect observed with TAROTs. 



Star 


P pu l B Known P B N max 


Ampl. 


New P B 


Full Amplitude 


Epoch [HJD-2450000] 


Phase shift 




[d] [d] 


[mag] 


[d] 


O-C [d] Mag. max. 


O-C Mag. max 


[P B \ 



Stars with known Blazhko effect 



RS 


Boo 


0.3773 


533 1 


88 


1.20 


532.481 





.012 


0.14 


3519.574 


3551.523 


-0.06 


ST 


Boo 


0.6222 


284 1 


58 


1.35 


284.090 





.032 


0.40 


1726.102 


1748.828 


-0.08 


AH 


Cam 


0.3687 


11. 2 4 


88 


1.06 


10.829 


0. 


.031 


0.29 


3632.301 


3635.117 


-0.26 


TT 


Cnc 


0.5634 


89 1 


82 


0.95 


89.02 


0. 


.020 


0.28 


2487.945 


2506.641 


-0.21 


Z 


CVn 


0.6539 


22. 75 5 


32 


1.00 


22.98 


0. 


.021 


0.32 


2745.629 


2763.781 


-0.79 


Z 


CVn 


0.6539 


22. 75 5 


32 


1.00 


22.98 





.021 


0.32 


2745.629 


2763.781 


-0.79 


ss 


CVn 


0.4785 


97 3 


39 


1.25 


93.72 





.028 


0.42 


3606.699 


3620.758 


-0.15 


RV 


Cet 


0.6234 


112 2 


23 




111.732 





.034 


>0.20 


0808.660 






RX 


Cet 


0.5737 


256 2 


15 




273: 


o 


061 




31 4^ Q99 






RX 


Col 


0.5937 


130 2 


20 


0.85 


134.77 


o 


.032 


0.29 


4188.273 


4228.703 


-0.30 


RY 


Col 


0.4788 






1 20 




n 
u 


022 


u.oo 


QflQS ^OA 




-U. 1 / 


VW 


Dor 


0.5706 


oa n2 
z4.y 


-1 1 

Tl 


1 on 


ot: no 

zo.yy 


n 




n oa 


4UyO.ZOZ 




n p.a 


RW 


Dra 


0.4429 


41. 61 1 


74 


1.30 


41.42 





.041 


0.25 


2371.344 


2374.656 


-0.08 


BD 


Dra 


0.5890 


24 3 


136 


1.25 


24.11 





.042 


0.40 


3605.965 


3607.894 


-0.08 


XY 


Eri 


0.5542 


50 2 


35 


1.35 


50.2 





.063 


0.75 


4143.031 


4162.109 


-0.38 


RX 


For 


0.5973 


31.8 2 


35 


1.20 


31.790 





.041 


0.24 


4080.688 


4087.680 


-0.22 


RR 


Gem 


0.3973 


7.230 6 


132 


1.25 


7.209 





.004 


0.11 


1936.137 


1939.957 


-0.53 


SZ 


Hya 


0.5372 


26.3 2 


85 


1.35 


26.233 





.053 


0.58 


2054.242 


2054.898 


-0.02 


IK 


Hya 


0.6503 


72 2 


9 


0.80 


67.5: 





.063 


>0.45 


4319.551 






RZ 


Lyr 


0.5112 


116.7 1 


101 


1.25 


120.190 





.017 


0.44 


0072.211 


0152.137 


-0.66 


SS 


Oct 


0.6219 


145 2 


28 


1.00 


144.93 





.007 


0.18 


4162.660 


4172.805 


-0.07 


uv 


Oct 


0.5425 


145 2 


30 


1.10 


146.99 





.012 


0.39 


5548.223 


5609.961 


-0.42 


RV 


UMa 


0.4681 


90.63 7 


87 


1.30 


90.170 





.005 


0.36 


0955.316 


0975.152 


-0.22 


AF 


Vel 


0.5274 


59 2 


30 


1.15 


58.548 





.017 


0.30 


4181.438 


4212.176 


-0.52 
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TY 


Aps 


0.5017 




25 


1.15 


109.13 





.016 


0.28 


4340.516 


4350.883 


-0.09 


BS 


Aps 


0.5826 




24 


0.80 


40.93 





.021 


0.24 


3166.379 


3182.137 


-0.38 


U 


Cae 


0.4198 




33 


1.35 


22.80 





.007 


0.25 


4459.238 


4473.602 


-0.63 


BI 


Cen 


0.4532 




35 


1.25 


79.68 





.021 


0.30 


4663.672 


4695.944 


-0.40 


uu 


Hya 


0.5239 




30 


1.20 


39.89 





.028 


0.26 


2367.930 


2401.035 


-0.83 


DD 


Hya 


0.5018 




25 


1.30 


34.59 





.016 


0.33 


2829.426 


2847.238 


-0.52 


RY 


Oct 


0.5635 




18 


1.25 


216.45 





.019 




4384.359 






ST 


Vir 


0.4108 




28 


1.30 


25.58 





.008 


0.10 


3484.652 


3495.016 


-0.40 


SV 


Vol 


0.6099 




31 


0.90 


85.47 





.009 


0.18 


4230.387 


4282.523 


-0.61 



Note . — Referenc es: HKholopov etaTI IU985I) : IWils fc Sodon 12005T ), TWils et al.l J2006I) , ILe Borgne et alj 1120071) , IKanvd 
(1966T) . " Uurcsik et all H2005aT) , ' IHurta et all 12003) 



tra, without having the possibility to determine individ- 
ual T max and V max . 

We analysed the OGL E stars considering th e different 
subclasses introduced by Collin ge et al.l J2006) on the ba- 
sis of the appearance of their amplitude spectra: BL1 (a 
single peak close to the main pulsational frequency), BL2 
(an equidistant triplet), and BL2+? (equidistant triplet 
plus an additional peak) stars. The distributions of the 
Blazhko periods vs. the pulsational ones are not very 
different in the plots of each subclass. The same applies 
for the similar categories used in the MACHO classifica- 
tion. Figure [12] allows a direct comparison, although the 
all-sky Blazhko variables are much less numerous than 
those in the galactic bulge and in the Large Magellanic 
Cloud. Blazhko periods longer than 400-500 d are not 
very common (see inserted boxes) . The lack of such long 
periods in the all-sky survey could be an observational 
bias due to the impossibility of following an isolated star 
along many years. 

The range of pulsational periods is quite similar and the 
very few Blazhko variables below 0.40 d and above 0.67 d 



reflects the lesser number of RRab stars having these pe- 
riods. The long Blazhko periods for stars with P < 0.40 d 
(RS Boo and of a couple of stars in the galactic bulge) 
are quite rare. Note how the limited all-sky sample en- 
hances an apparent linear increase of the Blazhko period 
from 0.40 to 0.55 d, accompanied by a rather flat dis- 
tribution starting from 0.47 d toward longest pulsational 
periods. This peculiarity is less discernible in the OGLE 
sample, but disappears in the MACHO one. We are not 
very confident of the reality of such structures. Indeed, 
the histograms of the Blazhko periods did not reveal any 
significant peak after the highest at 45 d in the MACHO 
sample and at 25 d in the OGLE one. 
We prefer to consider Fig. [TJ] as an evidence that the 
Blazhko period could be very different for the same 
pulsational period (vertical scatter) and that the same 
Blazhko period could be observed in a large range of 
pulsational periods (horizontal scatter). In such a con- 
text, the Blazhko effect of AH Cam fSect. | 3.1l) is very 
similar to that of VI 127 Aql (jChadid et alj|2010j) for the 
amplitudes of the magnitudes and phases of the maxima, 
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Fig. 10. — Stars with the brighest Vmax (ordinate) corresponding to the maximum negative O-C (abscissa). The closed curves are run 
in the anticlockwise direction. 



for the shape of the Blazhko cycle in the O-C vs. magni- 
tudes plot. The pulsational periods are very similar too. 
But the two Blazhko periods are very different, although 
both are short (Pb=10.829 d for AH Cam, P B = 26.88 d 
for V1127 Aql). 

The analysis of the LMC s tars in the MACH O sample 
allows another comparison. lAlcock et al.l (|2003[ ) report a 
higher incidence (74%) for frequency patterns in which 
the modulation amplitude on the higher frequency side- 
peaks is larger than its lo wer frequency coun terpart. The 
mathematical formalism ([Benko et al.ll201ll ) tells us that 
larger modulation on the higher frequency side indicates 
anticlockwise direction in the T max vs. O-C plots. As 
in the case of our galactic sample, this direction is the 
preferred one for LMC RRab stars. 

6. CONCLUSIONS 

The GEOS database and the TAROT survey provided 
a treasure of information on the Blazhko effect. TAROT 
light curves revealed that the non-repetitivity of the 
shoulder on the rising branch is a common feature in both 
Blazhko and non-Blazhko stars. This immediately trans- 
lates into a model constraint, since the physical cause 
creating the shoulder must work in both types of stars 
and the Blazhko effect is not damped by this cause. The 
hydrodynamical models should take into account several 
features (the bump at minimum, the shoulder on the ris- 
ing branch, and also the double maximum observed in 
RRc stars) that accurate photometry is putting in ev- 
idence on non-Blazhko stars, too. In this context, the 
new models expla ining the period doubling bifurcation 
(|Szab6 et al.ll2010l) should be used to test these features 
also. 

The Blazhko effect appears to be multifacetcd. We ob- 



serve very different shapes of the O-C vs. V max plots and 
the described paths (circular, elongated, almost linear, 
with loops and hooks, ...) can be run in a clockwise or 
anticlockwise direction also in stars with similar periods, 
both pulsational and Blazhko. Our observational results 
confirm the variety of shapes predicted by a mathemat- 
ical formalism invo lving amplitude and phase variations 
(jBenko et al.H20lTt ). As a further feedback to theoreti- 
cal works, the physical cause of the Blazhko mechanism 
must be able to explain both such a large variety and the 
preferred anticlockwise direction. The very-long photo- 
metric monitoring also suggests that the Blazhko effect 
shows cycle-to-cycle chang es (also in the case of RR Lyr 
itself; IPreston et al.lll965l Le Borgne et al. in prepa- 
ration) and this further variability has still to be im- 
plemented in the formalism. We could say that perfect 
regularity is the exception (i.e., BI Cen) rather than the 
rule. 

We emphasize the occurrence of different Blazhko pe- 
riods at the same pulsational periods and of the same 
Blazhko period at very different pulsational periods. We 
could be tempted to consider that the oblique-pulsator 
model combined with the randomly distributed line-of- 
sights and (very) different rotational periods could pro- 
duce the variety of the observed O-C vs. V max plots. 
However, the non-detection of relevant magnetic fields 
and the large cycle-to-cycle variations in the Blazhko ef- 
fect pose insoluble problems to this explanation. 
Since RRab stars are located in a narrow strip on the 
horizontal branch, their physical conditions are not ex- 
pected to vary too much. iLe Borgne et al.l (|2007f) showed 
how both redward and blueward shifts are observed in 
the evolutionary paths of RRab stars and perhaps the 
evolutionary stage (i.e., contraction or expansion) could 
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Fig. 11. — Same as Fig. 1101 but the closed curves are run in the clockwise direction. 



introduce relevant differences in the Blazhko mechanism. 
In the case of Z CVn we could verify that the varia- 
tion of the pulsational period resulted in changes in the 
Blazhko effect, both in terms of period and amplitude, 
and that they are anticorrelated. Since the variations 
of the pulsational period reflect changes in the star's 
radius, the Blazhko effect seems very sensitive to this 
stellar parameter. In this context, the unique case of 
correlated changes between the Blazhko and pulsational 
periods (i.e., XZ Dra), the few stars with P < 0.40 d 
showing a very-long Blazhko period and those showing 
complicated Blazhko plots can be considered promising 
laboratories in which to test future improvements in the 
theoretical models of the Blazhko effect. 
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and in the MACHO LMC survey (bottom panel). 
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APPENDIX 



THE INVENTORY OF GALACTIC BLAZHKO STARS 

The inventory of Blazhko stars locate d in the gala ctic field has been pro posed and revi sed several ti mes. The first com - 
prehensive list was published bv lSzeidll (119761 119881 ) and then adapted bv lSmithl (|1995f ) for his bo ok. iWils et al.l (l2006h 
discovered new Blaz hko RR Lyr stars by analyzing the Northern Sky Variability Survey (NS VS, iWozniak et al.l[2 004). 
Wils fc Sodorl (120051) did the ' same by analyzing the All Sky Au tomated Survey (AS A S, iPoimanski fc Macieiewskil 



2004). Sodor et al.l (|2005f ) corrected some un certain results an d Uurcsik et al.l (|2005bD added some complementary 



stars, partially overlapping those propo sed by IWils et al.l (|2006f ) . The Blazhko effect was discovered on many RRab 
stars located in the Kepler ficld-of-view (|Benk6 et al.ll2010l) . Table [2] lists the Blazhko RRab stars complem enting Ta- 
bleQ] New entries a re expected from the KBS II, e.g., V397 Her and ASAS 212034+1837.2 (jSodor et al.ll2012l), the Sloan 
Digit al Sky Survey (jSiiveges et al.l 127)121 ) . and from other specific projects, as SATINCP^I and FARO (|Le Borgne et al.1 
l2010h . 

Table[2]does not include the Blaz hko stars in the galac tic bulge discovered in the OG LE project, for wh ic h we refer the 
reader to the classification pape r (|Collinge et al.1120061 ). Moreover, the stars listed by IWils et al.1 (|2006l ). [Jurcsik et "all 
(2009b3), and lBenko et al.1 (|201Q[ ) with an uncertain Blazhko period are not considered as well. This is justified by the 



http: / /stiftung-astronomic.de/satino.htm 



12 



Le Borgne et al. 



59: d period sugges ted for RZ Lyr while the tr ue period is 120 d (lJurcsik et al.ll2012ft . Other cases are co ntroversial: 
FM Per (20 d after IWils et al.l l2006l 122 d a fterlLee fe Schmidtl |200l |), AR Ser (63 d afte r lWils et al.l l2006L 110 d after 
iLee fe Schmidtil2001l ). and X Ret (45 d after lSmithlll995l 161 d after IWils fe S6dorll2005ft . These uncertainties are not 
sur prising since sometim e s it is very difficult to have a clear evidence: t he Blazhko effect of DM Cyg was announced 
by ILvsova fe Firmanvukl (|1980n . then questioned by ISodor et al.l (|2005D . and finally confirmed by new observations 
(|Jurcsik et al.ll2009aj ) 

Except for the stars too faint to be measured with 25-cm telescopes, we are planning to complete their monitoring in 
the next years. In such a way we could extend the TAROT homogenous sample of well-studied Blazhko stars. 



Blazhko effect in RR Lyr stars 



TABLE 2 

RRab stars showing a known Blazhko effect and not discussed in this paper. 
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Not e. — See I Collinge et alJ 112009) for the OGLE s tars i n the galactic bulge. Ref erences: (1 ) ISmitbl 
(1995 1: (21 IWils fc S odor (2001 ) (3) ISodor et al.l (t2005T ). (41 IJurcsik et al.) l!2005bT) . ( 51 IWils et all (20061. 
(61 I Jurcsik et al.l H2009bl) , (7) ISzabo et al.l H20091) . (8) ILee fc Schmidtl lf200ir i~(9l IBenko et aTT j20icli . 
(10) Guescnbcreer et al. (2012). New names are provided by the General Cataloque_of Variable Stars for the 
newly discovered RRab Blazhko variables V1124 Her^NSV 8170 (IJurcsik et al. 2005b), PQ Lup^NSV 7330, 
DZ Oct^NSV 4350, CS P he=NSV 420, an d AL Pic^NSV 1700 {Wils & Sodor 2005), V1280 Ori=NSV 2724, 
and NS UMa^NSV 4034 dWils et al.| [2006D 



